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' The Schrodinger equation for quantum dot lattices with non-cubic, non-Bravais lattices built 

up from elliptical dots is investigated. The Coulomb interaction between the dots is considered in 
dipole approximation. Then only the center of mass (cm.) coordinates of different dots couple with 
each other. This cm. subsystem can be solved exactly and provides magneto- phonon like collective 
excitations. The inter-dot interaction is involved only through a single interaction parameter. The 
relative coordinates of individual dots form decoupled subsystems giving rise to intra-dot excitation 
modes. As an example, the latter are calculated exactly for two-electron dots. 

Emphasis is layed on qualitative effects like: i) Influence of the magnetic field on the lattice instability 
due to inter-dot interaction, ii) Closing of the gap between the lower and the upper cm. mode at 
' B=0 for elliptical dots due to dot interaction, and iii) Kinks in the intra dot excitation energies 

i-^ , (versus magnetic field) due to change of ground state angular momentum. It is shown that for 

lyj ' obtaining striking qualitative effects one should go beyond simple cubic lattices with circular dots. 

\ In particular, for observing effects of electron- electron interaction between the dots in FIR spectra 

(breaking Kohn's Theorem) one has to consider dot lattices with at least two dot species with 
different confinement tensors. 
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PACS: 73.20.D (Quantum dots), 73.20.Mf (Collective Excitations) 



I. INTRODUCTION 



Quantum dots have been in the focus of intensive research already for at least a decade which lead to a countless 
' number of publications [] (for a recent book see Ref. |l|). Although almost all experiments are performed at dot lattices, 
I> ' in the vast majority of theoretical investigations the interaction between dots is neglected. This is for the following 
' reasons: i) Because the confinement frequency is a parameter, which is mainly extracted from optical properties, it 
is difficult to tell the influence of dot interaction apart from the intrinsic single- dot value. (Possibilities to overcome 
this problem are discussed in the present work.) ii) The theory of Raman spectra, which can in principle monitor the 



dispersion (wave number dependence) of excitation energies as a direct consequence of interdot- interaction, is not yet 
advanced enough to extract the dispersion, iii) The lattice constant of dot arrays produced with current technologies 
, is so large (> 2000 A) that large electron numbers N per dot are necessary to obtain a seizable amount of shift. For 
^ ■ these N, however, reliable first principle calculations are not possible. With the advent of self- assembled dot arrays 

' the last item might change. 
^ r The scope of this paper is to investigate conditions, which lead to qualitative and observable effects of interdot- 
interaction on excitation spectra and the phase transition found in Ref. |[ Unlike in Ref. ||, a magnetic field B is 
explicitly taken into account and a microscopic theory is applied. Our approach is purely microscopic, i.e. we solve 
the Schrodinger equation of a model system exactly. Our model comprises the following approximations: i) The dot 
^ confinement is strictly parabolic in radial direction, but with anisotropic confinement frequencies LOi (i = 1,2) and 
|L* ' independent of N and B. ii) Overlap of wave functions between different dots is neglected (no hopping), iii) The 
. !^ ] Coulomb interaction of the electrons in different dots is treated in dipole approximation (second order in dot diameter 
over lattice constant). Our model is similar to that in Ref. ^ but allows more complicated dots and lattice structures. 
?H Besides we calculate also the intra dot excitations (apart from the collective center- of- mass excitations) for = 2 
explicitly and discuss the instabilitya in this microscopic model. Our results on the lateral dot dimer are compared 
with a former papercl, which uses a high magnetic field approach, in Sect. III. 

The plan of this paper is as follows. For further reference, we briefly summarize in Sect. 1 some relevant results for 
one single dot, or for dot lattices, where the distance between the dots is very large. This is important, because all 



^Therefore we will refer here only to papers which are directly connected to the scope of this work 



1 



exact solutions in the center- of- mass subsystem are traced back (by special transformations) to the solution of this 
one- electron Hamiltonian. This is analogous to ordinary molecular and lattice dynamics. After this, we consider 
a dot dimer, which mimics a lattice, where the dots are pairwise close to each other. This model can give an idea 
of the effects expected in dot lattices with a basis. Next we consider a rectangular, but primitive lattice in order to 
obtain the dispersion in the spectra. Finally, the intra-dot excitations of the Hamiltonian in the relative coordinates 
are calculated numerically. The paper ends with a summary. In the Appendix we give a short and elementary proof 
for the fact that the Generalized Kohn Theorem holds even for arbitrary arrays of identical non-circular quantum 
dots with Coulomb interaction (between the dots) in an homogeneous magnetic field. 



II. SINGLE DOT 



The Hamiltonian considered here reads (in atomic units ti = m 
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where m* is the effective mass (in units of the bare electron mass m), (5 the inverse dielectric constant of the 
background, and C a symmetric tensor. In case of a single dot, C is given by the confinement potential and we define 
C = n. It is always possible to find a coordinate system where VI12 = VL21 = and = Lof = m*Lo*'^ . We use the 
symmetric gauge A = ^ B x r throughout. The Zeeman term in H is disregarded at the moment. 
For iV = 1, the Hamiltonian 
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can be diagonalized exactly. Later on we will see that also the case of interacting dots can be traced back to the 
solution of type (|^). (Therefore, we kept the off diagonal elements of C in the results given below because the 
dynamical matrix, which also contributes to C, is generally non-diagonal and we want to use the same coordinate 
system for all q values.) After transforming thcj^perators r; and p; to creation- annihilation operators (see e.g. Ref. 
|l|) and using the procedure described by TsalliaJ, we obtain for the eigenvalues 
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'^o = 2(^11 + ^22) ; A = Cii-C22 (6) 

and Lo* ~ is the cyclotron frequency with the effective mass. (The results for the special case C12 = can also 
be found in Ref. |^.) The optical selection rules are the same as in the circular case, i.e., there are two possible types 
of excitations 

(An+ = ±1 and An_ = 0) or (An_ = ±1 and An+ = 0) (7) 

leading to the excitation energies /S.E = ujj^ and lo^ . It is easily seen that the form (^) reduces to the familiar formula 
in the circular case, where A = and C12 = 0. By inspection of (|^) we find that a soft mode uj-{B) = can only 
occur if Cii • C22 = C'^2- F'or a diagonal C this means that mm(Cii,C22) = 0. The last condition is of importance 
for interacting dots considered in the next Sections. 
In the limiting case B = we obtain from (|^) 
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u:±{B = 0) = \^ ^—^-'' ±^j ^-' +Cl, (8) 

We see that degeneracy uj+{B = 0) = lu^{B = 0) can only happen if C12 — and Cn = C22- For a diagonal 
confinement tensor with C12 = we obtain liJ+{B = 0) = max(cJi, 0^2) and ll!_(B = 0) = min(cji, tt'2)- As to be 
expected, we observe a gap between the two excitation curves uj+{B) and U!-{B) at i? = 0, if the two confinement 
frequencies do not agree. 

Alternatively we can introduce the quantum numbers 

(n+ - |n+ 

k ; = n+ - n_ (9) 

where k is the node number and turns in the circular limit into the angular momentum quantum number. 
For arbitrary N , the center of mass (cm.) R = can be separated H = Hem. + Hrei. with 
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where P = — iVr (see Appendix). As well known, Hem. does not contain the electron- electron- interaction. Hem. 
can be obtained from the one- electron Hamiltonian ^ by the substitution: B -> NB, C -> A^^C and H j^H. If 
we make the same substitution in the eigenvalues (|^), we obtain 

Ec.m.in+,n^) = EN=iin+,n^) 

i.e., the eigenvalues of the cm. Hamiltonian are independent of N. In other words, in H there are excitations, in 
which the pair correlation function is not changed, or classically speaking, where the charge distribution oscillates 
rigidly. Because FIR radiation (in the limit A 00) can excite only the cm. subspace, all we see in FIR spectra is 
the cm. modes. 

III. DOT DIMER 

We consider two identical elliptical dots centered at ai — (—a/2,0) and 02 = (+a/2,0). We expand the Coulomb 
interaction between electrons in different dots in a multi-pole series and restrict ourselves to the dipole approximation. 
By introduction of cm. and relative coordinates within each dot, the cm. coordinates and the relative coordinates 
decouple Q 

1,2 



H = Hem.iill,Il2)+J2 Hc,[{r}i''-'n (11) 



{r}a^ symbolizes {N — 1) relative coordinates in the a*'* dot. This means, we have 3 decoupled Hamiltonians: 
the cm. Hamiltonian and two Hamiltonians in the relative coordinates of either dot. This leads to two types of 
excitations: 

i) Collective excitations from Hem. which involve the cm. coordinates of both dots simultaneously. Because of 
the harmonic form (in the dipole approximation), there are exactly two modes per dot, thus a total of four. Each 
excitation can be classically visualized as vibrations of rigidly moving charge distributions of both dots. 

ii) Intra dot excitations which are doubly degenerate for two identical dots. Because i?Q,({r}i^ ^■') is not harmonic (it 
includes the exact Coulomb interaction between the electrons within each dot, which is not harmonic), this spectrum 
is very complex. It is the excitation spectrum of a single dot in a modified confinement potential where the cm. 
coordinate is fixed. The extra term in the modified confinement potential comes from the dipole contribution of the 
interdot Coulomb interaction. 

In this Section we consider only the cm. Hamiltonian and focus our attention to the the effects of ellipticity in the 



^ The tilde indicates that in this preliminary Hamiltonian a common gauge center for both dots is used. 
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dot confinement potential. The relative Hamiltonian for A'' = 2 is explicitly given in the last Section and solved for 
circular dots. For the elliptical confinement potential considered in this Section, the relative Hamiltonian cannot be 
solved easily, even if we restrict ourselves to = 2, because the elliptic confinement potential breaks the circular 
symmetry of the rest of the relative Hamiltonian. 



A. Center of Mass Hamiltonian of the Dimer 

The cm. Hamiltonian in the dipole approximation reads 
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where the small elongation Uq, is defined by Rq, = + and P = — zVr = —iVjj. The tensor C is 

Ca,a' = -PN T{aa,a') for a ^ a' 
where a^^a' — o,a ~ o-a' , and the dipole tensor is 

T(a) = ^ [iaoa-a^l] 



(12) 
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containing a dyad product (o) and the unit tensor I. As in the cm. system of a single dot, the explicit iV- 
dependence in (|l^) cancels in the eigenvalues. What is left is only the N- dependence in the dipole contribution of 
the dot interaction appearing in ( |l3| ) and (p^. This means, that the cm. spectrum of interacting dots is no longer 
independent of N . 

The term Oq, in the argument of the vector potential in (|2|) causes trouble in finding the eigenvalues. This shift is a 
consequence of the fact that we have to adopt a common gauge center for both dots (we chose the middle between 
both dots). This problem can be solved by applying the following unitary transformation 



i?c.™. = Hen. Q ; Q = n e"'* 



(16) 



In other words. Hem. agrees with Hem. except for the missing shift in the argument of the vector potential. 

The 4 modes inherent in Hem. are not yet explicitly known, because the 4 degrees of freedom are coupled. Decoupling 

into two oscillator problems of type (0) can be achieved by the following transformation: 
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and a is a vector pointing from one dot center to the other. Then T(a) has the following components 

Til ~ ; T22 = n ; T12 = T21 = 



(17) 
(18) 
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Now we assume that the principle axes of the confinement potentials are in x-y-direction. This means 



Qll = Luf ; ^22 = Uji ; ^12 = ^21=0 (22) 

The eigenvalues of can be obtained from (||) and (Q) with 

= ^(^? + ^2) ; A = (c^2 _ ^2) (23) 



and for iJ' ' with 



^o = li^l+^2} + lp; (ojl-^j) + 3p (24) 



where the interaction parameter is defined by 
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Observe that the dependence on N cancels, except that included in p (see discussion following (|10[)). 

It is important that the dot interaction influences the result only through a single parameter. This conclusion 
agrees with the semi- phenomenological theory in Ref. ^. 

In all our figures we express frequencies in units of the average confinement frequency luq = ^{uji +1^2), and A and p 
in units luq. Then, all systems can be characterized by the two parameters: uji/uj2 and p. In other words, all systems 
having the LO1/UJ2 ratio indicated in the figures are represented by the family of curves with the p values shown. The 
only exception we made is the cyclotron frequency lo*. ll!*/luo would be a good parameter in this sense, but we chose 
to use the magnetic field in Tesla instead for better physical intuition. The conversion between both scales is given 
by uj*[a.u.*] = B[Tesla] or, w*[wo] = §0^§i\^^B [Tesla] In this paper we used luq = 0.2 a.u.* = 2.53 meV 

and TO* of GaAs. (We want to stress that this choice effects only the magnetic field scale and not the qualitative 
features of the figures.) For easy comparison with experimental parameters we add the definition of effective atomic 
units {a.u.*) in GaAs (m* = 0.067, (3 = 1/12) for the energy: 1 a.u.* = 4.65 • lO"* double Rydberg = 12.64 meV , 
and for the length: 1 a.u.* = 1.791 • 10^ Bohr radii = 0.9477 • 10^ A. 

Because U'"'"-' agrees with the total cm. R = ^(Ri + R2) of the system, is the total cm. Hamiltonian. 

For B = 0, the eigenmodes can be visualized by classical oscillations. The two eigenmodes of are (rigid) in- 

phase oscillations of the dots in x and y direction, respectively. Because of the Kohn theorem (see Appendix), the 
independence of iJ*^"*"^ on the Coulomb interaction does not only hold in the dipole approximation, but it is rigorous. 
This shows also that the dipole approximation is consistent with the Kohn theorem, which is not guaranteed for single 
particle approaches. Because FIR radiation excites (in the dipole approximation) only the cm. modes, it is only the 
p-independent eigenmodes of which are seen in FIR absorption experiments. 

This statement is in contradiction to Ref. |^. They performed numerical diagonalizations for a lateral pair of circular 
dots confining the set of basis functions to the lowest Landau level and considering parallel spin configurations only. 
This is justified in the limit of high magnetic fields. They found a splitting of the two dipole allowed modes at i? = 
due to dot interaction and some anti-crossing structures in the upper mode, whereas the lower mode is always close 
to the single particle mode. This fact is already a strong indication that the missing higher Landau levels cause both 
spurious effects. (Observe that the lifting of the degeneracy at i? = in the dipole allowed excitations in Fig.l is due 
to the ellipticity of the intrinsic confinement and not due to dot interaction.) 

The eigenvalues of H^~^ do depend on p because the dots oscillate (rigidly) in its two eigenmodes in opposite phase, 
one mode in x and one mode in y direction. This leads to a change in the Coulomb energy. The two eigenmodes of 
can also be described as a breathing mode (in x direction) and a shear mode (in y direction). 
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Dot dimer, (Oj/o)j=3/2 
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FIG. 1. Excitation energies AE^"^ = oj^'' for a dimer of elliptical dots as a function of the magnetic field for some discrete 
values of the interaction parameter p {p values in the inset axe in units ljq). The ratio of the oscillator frequencies in the 
direction of the capital axes loi/lu-a is a) 3/2 and b) 2/3. The dipole allowed excitations AS^"*"^ of H^'^^ (thick full line) are not 
infiuenced by the dot interaction and therefore independent of p. 



B. Special Features of the Excitation Spectrum 

In Fig. la and lb, the four excitation frequencies of the dimer are shown with p as a parameter. 
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For p = 0, the two modes ^ agree with the two modes lo^±^ . In all symbols, the superscript sign refers to the 
system H<-+^ and H<^-^ (c .m. or relative coordinate), and the subscript sign discriminates the two modes of the same 
system. The two modes w^"^ are independent of p. There are two qualitatively different cases. (Consider that toi 
is the oscillator frequency parallel to the line, which connects the two dot centers, and lu2 is the oscillator frequency 
perpendicular to this line.) If uJi > lo^ (Fig- la), the gap between and wL"'' at _B = increases steadily with 
increasing p until, for a critical Pcr = ^\ (iu our numerical case: pcr['^o] ~ 16/25 ~ 0.64) the lower mode u}_ ^ becomes 
soft. This transition is independent of B. For o^i < ui the gap between cj^ ' and iJ'_ ' at i? = first decreases with 
increasing p until it vanishes for p ~ ^(t^i ^ "-^i) (i^^ o^'' numerical case: p\lo'^ = 4/15 = 0.27). Afterwards, it increases 
until the lattice becomes soft at pcr = UJ2 our numerical case: pcr-[wg] = 36/25 = 1.44). The dependence of the 
two excitation energies w^"'' and lu^S'' on p for B = in the second case is shown in Fig. 2. Comparison of Fig.s 2a 
and 2b demonstrates that the dot architecture in Fig. 2a is much more sensitive to interdot interaction than that in 
Fig. 2b. Thus, if we want to observe or use the instability, this event happens in case 2a for for smaller p (or larger 
lattice constants) than in case 2b. Additionally, the assumption of non-overlapping dot wave functions (for a given 
lattice constant) is better fulfilled in case 2a than in case 2b. 

Dot dimer, co,/co =2/3, B=0 




FIG. 2. Excitation energies AE\_ ' — tj£ ' of the Hamiltonian H'- ' for B = and Lo^juji = 2/3 as a function of the 



interaction parameter p. 



For GaAs as a typical substance, (E5|) can be rewritten in more convenient units as 



2.26- 10^ N 

(a[l])3 {u;o[meV]y 



Obviously, we need large dots (large N , small ujq- which means large polarizability), and a small dot distance a for 
seizable interaction effect. On the other hand, the dot radius for iV = 1 is of order of the effective magnetic lengt 

^0 ~ (2ti'o)^ + (^*)^ ) , which reads for GaAs 

lo [A] = ^-^ — (27) 

[LUo[meV]f + 0.739 {B[Tesla]y 
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and we need small dots and high magnetic fields for small overlap. Consequently, a magnetic field helps avoiding 
overlap of the dots, although e.g. the critical p for soft modes is independent of B. The question is, if there exists 
a window between these two (partly) conflicting demands. For an order- of- maigntude estimate, let us consider 
GaAs with cjo as chosen above and the worst case N = \. Then (|2^ ) with a typical dcjg] =0.1 (which seems to be 
the minimum for any observable effect) provides a dot distance of a [A] — 327 and (E^ gives for £? = a radius of 
Iq\A\ = 150 and for B\Tesla] = 10 a radius of Iq[A\ — 80. Consequently, the constraint /o < a/2 for our model can be 
fulfilled. For obtaining larger interaction effects the parameters have to be optimized. 

The next question is what happens in mode softening physically? Firstly, it is the antisymmetric shear mode lo'i_ ' 
which has the lowest frequency and which becomes soft. If the interaction parameter is strong enough [p > Pcr), 
the decrease in interdot- Coulomb energy with increasing elongation of the dots becomes larger than the increase 
of confinement potential energy. Because in the harmonic model both energies depend quadratically on elongation, 
the dimer would be ionized, i.e. stripped of the electrons. Clearly, in this case we have to go beyond the dipole 
approximation for the interdot interaction and beyond the harmonic approximation for the confinement potential. 

In order to obtain a hand-waving picture of what happens, the confinement potential of the system for shear mode 
oscillations is supplemented by a 4*'* order term in the following way: Vconf. — 2iV [-^1^2 — A [/^] with {A > 0), 
and the Coulomb interaction in 4*'* order reads: Vmt — —pN + (SpN/a^) where p — 2N(3/a^ as above. Then, 
the stability condition reads 



Yl^^iu;i-p)U' + i^-2A)U'>0 (28) 



The condition for the existence of a bound state is that the [/^- term is positive: Zp/a? > 2 A. For a positive U^- term 
( p < UJ2 )j the equilibrium position is C/q = 0. If the U'^- term becomes negative ( p > w| ), the system finds a new 
equilibrium at a finite elongation 



C/„ = ±/jpiiL (29) 

This new ground state is doubly degenerate: Ui = {—a,+Uo) , U2 = {+a,—Uo) and Ui — {—a,—Uo) , U2 = 
i+a, +Uo) have the same energy. In short, at Pcr = w| there is an electronic phase transition to a polarised state, 
where the equilibrium position of the cm. is no more in the middle of the dots. At the end we want to stress that 
all these stability considerations are only valid if the confinement potential is not changed under elongation of the 
cm. of the dots. Secondly, it is not rigorous to include the fourth order terms after separation of cm. and relative 
coordinates, because in fourth order these two coordinates do not decouple exactly. 

IV. DOT LATTICE 

We consider a periodic lattice of equal quantum dots at lattice sites rI^L = Rn''' + cia- The vectors R^"' form a 
Bravais lattice and runs over all sites within an unit cell. In developing a theory for these lattices we have to 
repeat all steps in Sect. II from ( pT| ) to (|l^) just by supplementing the index a by the index n for the unit cell. 



A. Center of Mass Hamiltonian of the Dot Lattice 

The cm. Hamiltonian in the dipole approximation then reads 
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where XJn.a — Rn.a — Rn^a is the elongation of the cm. at lattice site (n, a) and the force constant tensor C is 
defined in analogy to (O) and ([1^). 

The Hamiltonian (|30|) is a phonon Hamiltonian in an additional homogeneous magnetic field. The first stage of 
decoupling can be achieved by the usual phonon transformation 
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(32) 



where Nc is the number of unit cells and the transformed coordinates have the following properties U_q_Q, — U* ^ = 
Q, and P-q,a = Pq a- Thc Hamiltonian in the new coordinates is a sum of iVc subsystems of dimension 2x number 
of dots per unit cell: Hem. — X]q -^qi where 



^ TV 1 ^ 2m* 



Pq,„ + -A(U;j 



Pq,„ + -A(u; J 



^ ^ Uq ^ • Cq;Q,_Q' • Uq 



The dynamical matrix is defined by 



(33) 



(34) 



and it is hermitean Ca^^a' ,a — Cq.^, = C_q;Q_Q'. 

Next we want to recover the limiting case considered in Sect. III. If the dots in a given unit cell are far away from 
those in neighboring cells, then in (34) only the term with r'i^'' — contributes, C does not depend on q, consequently 
the index q is redundant, and (|3^) agrees with (|l2|). 

Our preliminary result ( p3| ) is not yet diagonal in a. a' . In some special cases (see e.g. two identical dots per unit cell 
considered in Sect. Ill) this can be achieved by an unitary transformation 



Uq^Q, ^ ^ Qq;a,Q' ' Uq , 



?q;Q',Q Qc^■a,a' 



(35) 



under which the one- particle term in (33) is invariant and the transformed interaction term 

-Z E ^q Q ■ Cq-a,a' ' Uq,Q' with Cq-^a.a' ~ ^ Q a]ai ^q:Ql, 



,a2 ^CK2,CK 



(36) 



can be made diagonal Cq;^ ^' = C^-a Sa,a' by a proper choice of Qa,a'- Now, (|3^) reads Hq = ^q,ai where 
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q,ct 1 ^ A(Uq_Q,) 



TT* . P . TT 
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(37) 



The eigenvalues of ( p7| ) can be obtained from those of (g) because corresponding quantities have the same commuta- 
tion rules. Such an unitary transformation does not exist, e.g., for two different dots per cell. Then (^3|) has to be 
solved directly using the method described in Ref. 



B. Dynamical Matrix for Bravais lattices 

From now on we consider Bravais lattices what means that we can forget the indices a in the first part of this 
Section. Then the dynamical matrix 

Cq = ^l + pN J2 (l - e^''''"" ') T (Ri°^) (38) 
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is real and symmetric, but generally not diagonal, even if n is diagonal. A very important conclusion is apparent 
in (^). In the limit q — > 0, the inter- dot interaction (represented by (5) has no influence on Cq and therefore on 
the spectrum. This means, that the excitation spectrum observed by FIR spectroscopy is not influenced by inter- dot 
interaction and agrees with the one- electron result (as in the single dot). This statement is rigorous for parabolic 
confinement (see Appendix). It can also be understood intuitively, because a g = 0- excitation is connected with 
homogeneous in- phase elongations of the dots which do not change the distance between the electrons. We want 
to mention that this conclusion seems to be in contradiction with the experimental work in Ref. |ll|. They found 
a splitting of the upper and lower excitation branch at _B = and g = for circular dots in a rectangular lattice, 
which they interpreted within a phenomenological model of interacting dipoles as a consequence of lattice interaction. 
However, they use mesoscopic dots with a diameter of 370000A and lattice periods of 400000 and 8OOOOOA. These 
dots are clearly beyond our microscopic quantum mechanical model, which rests on a parabolic confinement. 

For the rectangular lattices considered in our numerical examples we define R^*^) — NiOiei + N2a2G2 aud q = 
qi^ei + 52^62 with the lattice constants oi and 02 and integers A^i and A^2 characterizing the lattice sites. The 



components of q vary in the Brillouin zone (BZ) in the range [—1/2 
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T(iVi,7V2) = 
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1/2]. The dipole tensor ( [l5| ) reads 
3iViiV2aia2 



{2Nlal - Nlal) 



(39) 



Although for all figures the exact dynamical matrix is used, it is useful to consider the results with nearest neighbor 
(n.n.) lattice sums in (Bq) separately. This provides simple formulas for order- of- magnitude estimates. 



Cii = ul H 
C22 = ^2 H 

C12 = rii2 



2pi [1 — cos(27rgi)] — p2 [1 — cos{2'Kq2) 
2p2 [1 - cos{2Trq2)] - pi [1 - cos{2nqi) 



(40) 



where we introduced the interaction parameters pi 
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(for a,=2a^ in the middle of the irreducible Brillouin zone) 



e- - 



-e o 



e — e s„ 

G--G 

G---0 s,2 

Ofor nearest neighbour sums 



G 0. 



^-eo — o- 



-0 — 




10000 



FIG. 3. Convergency of the lattice sums Sik as defined in ( |4l[ ) with increasing number of cubic shells Nmax for q in the 
middle of the Brillouin zone. 



The convergence of the lattice sums Sik m the dynamical matrix is shown Fig. 3. Sik is defined by 

Cik — ^ik + P2 Sik 



(41) 
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and depends only on q and the ratio axjai. Apart from the ofT-diagonal elements, which vanish in n.n. approximation, 
the error of the n.n. approximation is less than 30%. 

C. Special Features of the Magneto— Phonon Spectrum 

Fig.s 4-6 show the magneto- phonon ^ spectrum for circular dots on a rectangular lattice with a\ = la^- Because 
the two interaction parameters have a fixed ratio, it suffices to use one of them for characterizing the interaction 
strength. We chose the larger one p2 = V- For 5 = and isolated dots {p — 0), the two excitation modes are 
degenerate. If we tune up the interaction strength represented by p, a q dependent splitting appears (see Fig. 4). 



B=0 



p[«„ ]=o 

0.25 - 

0.5 






'\ 


/ / ' 


\ 

\\ 




^y j/ 







I ' ' ^^-^ 1 

(0,0) (1/2,0) (1/2,1/2) (0,1/2) (0,0) 

(qi.qs) 



FIG. 4. Magneto- phonon dispersion at i? = for several interaction parameter values and q on symmetry lines of the 
Brillouin zone, (p = p2 values in the inset are in units tjQ.) Thick (lower) and thin (upper) lines indicate Ai?- and Ai?+, 
respectively. 



This splitting is a manifestation of the dot interaction. For a certain critical pcr the lower mode becomes soft. 
This feature will be discussed below. There are points in the BZ, however, where the degeneracy for finite p remains. 
These points will be investigated now. We demonstrated in Sect. II after formula (|8|) that necessary for degeneracy is 
Ci2 ~ 0, i.e., the dynamical matrix must be diagonal. Then the points with degeneracy are defined by the condition 
Cii = C22- As seen in (^Sf), for circular dots ui = L02 = this happens in the center of the BZ q = 0. The next 
question to be discussed is if there are other points with degeneracy. The first condition C12 = 0, is fulfilled for 
all points on the surface of the BZ. The second condition must be investigated for special cases. We find, that for 
quadratic lattices ai — 02 with circular dots uji = 102 both modes are degenerate at the point q = (1/2, 1/2). In the 
case shown in Fig. 4 this point is somewhere between (1/2, 1/2) and (1/2, 0). 

In n.n. approximation (^0|), however, this equation is even fulfilled on full curves in the BZ defined by 
Pi [1 — cos (2Trqi)] = P2 [1 — cos(27rg2)]- In a cubic lattice, this is the straight lines q2 = iqi- The contributions 
beyond n.n.s remove the exact degeneracy on this curve in the interior of the BZ, but leave a kind of anti-crossing 



^ The term magneto- phonon is attributed to the fact that the there is no exchange and there are harmonic forces between 
the oscillating individuals. One could also call them magneto- plasmons, if one wants to emphasize that it is only electrons 
which oscillate, and no nuclei. 
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behavior of the two branches. 



An important parameter, which characterizes the influence of the dot interaction in circular dots, is the band width 
at B — 0, i.e. the maximum sphtting of the two branches due to dot interaction. (Remember that this sphtting 
vanishes for noninteracting circular dots.) Assume oi > a2. Then the largest splitting for circular dots in n.n. 
approximation appears at q = (0, 1/2) and has the amount 



W ~ maa;(A_E_|_ — AE^) = max{Lu^ — lu^) = 



For small dot interaction and in units luq, this is proportional to the interaction parameter — 3p2 



(42) 



We next discuss the appearance of soft modes. The question is, for which q, B and interaction parameter p this 
happens. The general condition for vanishing of the lowest mode is Cu ■ C22 — C12 (see Sect. II). In this condition 
the magnetic field does not appear. For circular dots and with the definition ([41]) this equation reads 



[ujI + P2 Sii] [ujI + P2 Sii] = pI SI2 



(43) 



After introducing a dimensionless critical interaction parameter Per = PiI'^'q^ we obtain a quadratic equation for 
which has the solution 



P., 



1 tr 
2IM 




(44) 



where det = S11S22 — 'S'12 s-nd tr = Sn + 822- For our numerical case ai = 2 02 and n.n. interaction for Sik the 
lowest mode becomes soft at q = (0,1/2) and the critical interaction parameter is Per = 1/2. Inclusion of lattice 
contributions beyond n.n. shifts Per to 0.7543. The most important result of this paragraph is that lattice softening is 
independent of B (see also Fig.s 5 and 6). The latter conclusion is exact within the range of validity of the Hamiltonian 
(|^) and no consequence of any subsequent approximation or specialization. 

P=Por 




(0,0) 



(1/2,0) 



(1/2,1/2) 



(0,1/2) 



(0,0) 



FIG. 5. The same as Fig. 4, but for p = pcr and several magnetic fields. 
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V. INTRA-DOT- EXCITATIONS FOR N=2 



Intra-dot excitations for circular dots in a cubic lattice and for N = 2 can be calculated easily. We define the 
relative coordinate r = r2 — ri, and assume that all dots are equivalent (also with respect to their environment). Then 
the indexes (n, a) can be chosen as (0, 0) and omitted. The relative Hamiltonian reads 

+ lr.D.r+A| (45) 
where p = — zVr and 

D = ln + ^To; To= Yl t(r(°L) (46) 

It is worth emphasizing that Hrei contains a contribution from neighboring dots, originating from the interdot Coulomb 
interaction. A trivial angular dependent part can only be decoupled from Hrei, or, the 2-dimensional Schrodinger 

equation can be traced back to an ordinary radial Schrodinger equation, if the term r • D • r has the same circular 
symmetry as the intra-dot Coulomb term /3/(2r). Therefore we confine ourselves to circular dots on a cubic lattice, 
and we have 

where the simple result is in n.n. approximation. Using the interaction parameter p = 2N (with N = 2) defined 
above, we obtain 

D = ^(a;2 + 2p)I (48) 



Hrei 



2m* 



P+^A(r) 
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In this way, dot interaction defines an effective confinement frequency uJq 



2p. This means that the cm. 



excitations have to be calculated (or interpreted) with another confinement potential then the relative excitations. In 
our figures we present results for ti^o.e// = 0.2 a.u.* , which agrees with the bare confinement potential used in Sect. IV 
and the mean value in Sect. III. Because our results are presented in units of loq, they depend on too only weakly 
through the differing influence of electron-electron interaction. For the absolute values, however, the influence of the 
dot interaction can be tremendous. 

In the relative motion there is a coupling between orbital and spin parts through the Pauli principle. For N = 2 
and a circular effective confinement, Pauli principle demands that orbital states with even and odd relative angular 
momentum rui must be combined with singlet and triplet spin states, respectively ( see e.g. Ref. For the cm. 
motion there is no interrelation between orbital and spin part because the cm. coordinate is fully symmetric with 
respect to particle exchange. Consequently, any cm. wave function can be combined with a given spin eigen function. 
The only spin dependent term in the total energy considered here is the Zeeman term, which reads in our units 



— = 0.9134- 10"^ 

LOO 



B[Tesla] Ms 
LOola.u.*] 2 



(49) 



where we used Qs — —0.44 for the gyro-magnetic factor of GaAs from Ref. ^. The total spin quantum number is 
Ms = for the singlet state and 0, ±1 for the triplet state. One of the most interesting points in quantum dot physics 
is that the total orbital angular momentum of the ground state depends on the magnetic field (see e.g. Ref.s ^, |^). 
This feature is a consequence of electron- electron interaction. For our parameter values, the relative orbital angular 
momentum of the ground state changes from to -1, from -1 to -2, and from -2 to -3 at B = 1.250, 4.018, 
and 5.005 Tesla. This corresponds to a sequence Ms=0,+1,0,+1 for the spin quantum number. Figures 7a-c show 
the excitation frequencies for three B-values lying within the first three regions, is the relative orbital angular 
momentum of the final state. All excitations are included irrespective of selection rules. For dipole transitions only 
two of them would remain (the lowest excitation with riij = to^ ± 1). For B ^ 0, the lowest excitation energy (in 
units luq) for noninteracting electrons would be 1. As seen in in Fig. 7a, electron- electron interaction decreases this 
value by at least a factor of 1/2. The same holds qualitatively for finite B. This is connected to the fact, that the 
ground state depends on B. Let's consider an example. For B — 1.250 Tesla the ground state switches from mi = 
to rrii = —1. This implies that for B approaching this transition field from below, the excitation energy for dipole 
allowed transition from mi = Otom/ = — 1 converges to 0. In other words, there is a level crossing at the the 
transition field. Therefore, very small transition energies and switching of the ground state are connected. 



4.0 



intra-dot excitation energy, B=0 

(final state m,, initial state m.=0) 



3.0 - 



^ 2.0 

< 



1.0 



0.0 





m, 
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4.0 



3.0 



^ 2.0 
M 

<J 



intra-dot excitation energy, B=3 Tesia 

(final state m,, initial state m=-^) 



1.0 



0.0 



-4 -3 -2 -1 



4.0 



3.0 



^ 2.0 

< 



intra-dot excitation energy, B=4.5 TesIa 

(final state m,, initial state m =-2) 



1.0 



0.0 



FIG. 7. Intradot excitation energies (in units of the effective confinement frequency) for B=0 
TesIa (c). The corresponding relative orbital angular momenta of the ground state are mi=0 (a), 
angular momenta Ms=Q (a), +l(b), and 0(c). 



a), B=3 Tesla (b) and B=4.5 
-1(b), and -2(c) and the spin 



For a qualitative understanding, Figures 7a-c can be used together with Figures la, lb, 4, and 5 to investigate the 
relative position of collective and intra-dot excitations. The conclusion is that for small dot interaction (for p well 
below Pcr)i the lowest intra-dot excitation energies lie well below the lowest cm. excitations. Apart from using a 
different terminology, this conclusion agrees with the experimental findings in Ref. |l^. 

Fig.s 7b and c, which belong to finite -B, show the Zeeman splitting. All transition energies to final states with odd 
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m / are triplets because the corresponding spin state is a triplet state. The thin lines of a triplet belong to spin- flip 
transitions. 

intra-dot excitation energy 



(from ground states with m, to states with m,) 




0.0 1.0 2.0 3.0 4.0 5.0 



B [Tesia] 

FIG. 8. Intradot excitation energies (in units of the effective confinement frequency) as a function of the magnetic field. 
The small Zeeman splitting is neglected. The B- values, where the angular momenta of the ground state change, are indicated 
by vertical lines. The absolute value of the final state orbital momentum ruf of the curves at B=0 grows from bottom to top 
by 1 starting with 0. 



In Fig. 8 the B dependence of the lowest excitation energies is shown. It is clearly seen that the curves exhibit a kink 
at those B- values, where the ground state configuration changes. The size of the kink decreases with increasing B. 
If this kink could be resolved experimentally (e.g. by electronic Raman spectroscopy), it would be a direct indication 
for the change of the ground state configuration, and thus an experimentally observable consequence of electron- 
electron interaction. 
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VI. SUMMARY 



We solved the Schrodinger equation for a lattice of identical parabolic (but not necessarily circular) quantum dots 
with Coulomb interaction (in dipole approximation) between the dots. We provide an overview over the state of art 
of these systems which includes the results of former publications. References can be found in the text. 

• Similar to single dots, the center of mass coordinates of all dots can be separated from the relative coordinates. 
Only the cm. coordinates of different dots are coupled to each other. The relative coordinates of different dots 
are neither coupled to each other nor to the cm. coordinates. 

• This gives rise to two types of excitations: two collective cm. modes per dot and and a complex spectrum of 
intra-dot excitations. In periodic arrays only the collective cm. modes show dispersion. Intra-dot excitations 
are dispersion-less. 

• The cm. system can be solved exactly and analytically providing magneto- phonon excitations characterized 
by a certain wave number q within the Brillouin zone. For q = and one dot per unit cell, interdot interaction 
does not have any influence on the cm. excitations. 

• All dipole allowed excitations (seen in FIR experiments) are not influenced by the dot interaction. 

• Interdot interaction between two dots influences the spectrum through a single parameter p = 2Nf3/a^, where 
a is the distance between the dots, N the number of electrons per dot and /? the inverse background dielectric 
constant. 

• If p exceeds a certain critical value Pcr, the lowest cm. mode becomes soft leading to an instability. This 
transition is independent of the magnetic field. 

• For B=0 and and one circular dot per unit cell, the two cm. modes are not only degenerate in the middle of 
the Brillouin zone, but also at some points on the surface. If we use the n.n. approximation for the lattice sums 
in the dynamical matrix, degeneracy is maintained even on full curves in the Brillouin zone. 

• Intra-dot excitations have to be calculated from an effective confinement. In circular dots with a cubic envi- 
ronment in nearest neighbor approximation the effective confinement frequency reads uJq ^ = ojq + 2p. This 
effective confinement differs from that for the cm. motion. 

• For p well below pcr-, the lowest intra-dot excitations are much smaller than the lowest collective excitations. 

• The intra-dot excitation energies versus magnetic field exhibit kinks at those fields, where the angular momentum 
of the ground state changes. 

In the Appendix we prove a Kohn Theorem for dot arrays with Coulomb interaction between the dots without the 
dipole approximation. The individual confinement potentials can be arbitrarily arranged and can carry different 
electron numbers, but have to be described by identical confinement tensors. This means that for breaking Kohn's 
Theorem in dot arrays, we have to have at least two different confinement species. 



APPENDIX 

We are going to prove that for an arbitrary array ^ of identical parabolic quantum dot potentials in an homogeneous 
magnetic field: i) the total cm. degree of freedom can be separated from the rest, ii) the total cm. Hamiltonian is 
not influenced by Coulomb interaction, and iii) the eigenvalues of the total cm. Hamiltonian are independent of the 
electron number N in each dot ^ . 

The Hamiltonian H = -|- V consists of an one-particle term iJ = and the Coulomb interaction between all 



*The dot centers can be arranged arbitrarily. 

^The confinement tensors f2 of all dots must be equal. 

® The electron number in different dots can be different. 
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electrons V. The dot centers are located at R^"^ and the electron coordinates are denoted by = Ri°^ + ■ Then 
we have 



'^-^ I 2m* r- 



(50) 



First of all we shift the gauge center for each electron into the middle of the corresponding dot using an unitary 
transformation similar to (|l^). This transforms the shift 'Ra'^ in the argument of the vector potential away. Next we 
drop the index a in ( |50| ) so that the index 'i' runs over all electrons in all dots. Now we perform a transformation to 
new coordinates u,- 



u. = {Vn Uk) ; (ViV u,) = 



(51) 



where Qik is an unitary matrix. This implies 



(52) 



It is possible to choose for the first column Qki — The other columns need not be specified. Then Ui = 

(1/iV) Uj = U is the cm. of all elongations, or, the cm. of the electron coordinates with respect to the weighted 

center of the dot locations R^"^ = (1/^) J2a Rq'"' , where Na is the number of electrons in dot a. The corresponding 
canonical momentum pi — (l/z)Vai = P is the cm. momentum. Inserting our transformation into (^ ) provides 



^ I 2m* 



N 



A(ii,) 



Ui • C • U; 



(53) 



The term z = 1 in ( p3|) is the (separated) c.m. Hamiltonian 



2m* 



N 



A(U) 



2 



(54) 



which agrees with (10). Clearly, the Coulomb interaction F in is independent of the cm., and docs not contribute 
to Hem.- For the independence of the eigenvalues of N see the discussion following (p^). 



This prove, in particular the step from (|5C|) to (|5^), is not correct if the dot confinement tensor C depends on a 
(or in the changed notation). Therefore, all dots must have the same C, but can have different electron numbers 
Na- In other words, the total c.m. excitations in dot arrays, which are seen in FIR spectra, are not affected by the 
e- e- interaction, if and only if all confinement tensors C are equal. On the other hand, if we want to observe e- 
e- interaction in the FIR spectra and break Kohn's Theorem, we have to use dot lattices with at least two different 
confinement tensors. The simplest way to implement this is using a lattice with two non-circular dots per cell, which 
are equal in shape, but rotated relative to each other by 90 degrees. 
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